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When pairing correlations in a quasi two dimensional electron system induce a pseudogap in the single particle 
density of states, the specific heat must also contain a sizeable pair contribution. The theoretically calculated 
specific heat for such a system is compared to the experimental results of Loram and his collaborators for un- 
derdoped YBa2CusO^+ x and La,2- x Sr x Cu04, samples. The size and doping dependence of the extracted 
pseudogap energy scale for both materials is comparable to the values obtained from a variety of other experi- 
ments. 
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Introduction. The gradual, but substantial loss of low en- 
ergy single particle states in the normal state of underdoped 
cuprate samples is by now documented by a large num- 
ber of experiments [|J. The onset of the pseudogap regime 
is captured most clearly by spectroscopic probes that cou- 
ple predominantly to the single particle excitations, such as 
tunneling \^ and angle-resolved photoemission spectroscopy 
(ARPES) [[|]. All theoretical expectations, however, point 
towards some sort of many-particle phenomenon behind the 
pseudogap effect, and imply observable changes in collec- 
tive electronic properties [Q], Traditional probes of collec- 
tive degrees of freedom, such as inelastic neutron scatter- 
ing Raman scattering [g], optical conductivity [Jj]], were 
used extensively on pseudogapped samples, while new exper- 
iments, such as higher order tunneling spectroscopy [0], were 
also proposed. Despite such remarkable accumulation of high 
quality experimental data a consensus has yet to emerge re- 
garding the origin of the collective effect causing the pseudo- 
gap- 

One of the earliest experimental indications of normal state 
anomalies in underdoped cuprates came not from spectro- 
scopic but from thermodynamic measurements. Loram and 
his collaborators found [||] that the coefficient of the electronic 
heat capacity 761(1") = C e \(T)/T of underdoped samples is 
no longer constant in temperature, as it is for the optimal and 
overdoped crystals, and as one would expect from a normal 
Fermi liquid. Instead, 7 i(7 1 ) shows a broad maximum at a 
crossover temperature T* which is doping dependent and typ- 
ically much higher than the superconducting transition tem- 
perature T c . As T is lowered below T*, j e \(T) decreases sig- 
nificantly before the temperature reaches T c . Despite the fact, 
that these specific heat measurements have been available for 
some time, there is no systematic study of this data within any 
theoretical framework of a proposed pseudogap scenario. The 
purpose of this paper is to provide such an analysis within a 
pairing fluctuation framework. 

There are several reasons why the results of these ther- 
modynamic measurements deserve a detailed theoretical 
analysis. First, there is ample experimental data on the 
most extensively studied classes of cuprate superconduc- 
tors, YBa2Cu30e+ x , La 2 - x Sr x Cu04 and recently 



Bi^Sr-iC aCu-iOz+x for a remarkably wide range of dop- 
ing. This makes it possible to single out and to concen- 
trate the theoretical study on the universal features mentioned 
above. Second, these thermodynamic measurements couple to 
all thermally excitable modes, irrespective of their single par- 
ticle, or collective character. In tunneling and photoemission 
the higher order excitations, such as pair excitations, may be 
detected only through their convoluted influence on the single 
particle spectrum. Thus, in contrast to any theoretical analy- 
sis of spectroscopic data, where one seems to have a choice 
in selecting the collective phenomena causing the pseudo- 
gap, the comparison between theory and thermodynamic data 
provides a consistency check, once the choice has been made. 

In this article we perform such a consistency check for a 
simple pair fluctuation model of the pseudogap regime, intro- 
duced by Vilk and Tremblay [|To|]. We opted for this formu- 
lation because it permits essentially analytical treatment, and 
therefore the analysis is not obscured the numerical difficul- 
ties ubiquitous in most other incarnations of the pair fluctua- 
tion scenario [|i"T}]. We believe however that the results we ob- 
tained have a range of relevance that goes beyond the present 
formulation and it is characteristic of a larger class of pah- 
fluctuation models for the pseudogap [Q. 

In our calculation we consider beside the usual single parti- 
cle contribution [p3[], the contribution of the fluctuating pairs. 
We show that (a) it has the same order of magnitude as the sin- 
gle particle contribution and, (b) it is essential for explaining 
the universal presence of a broad hump in the specific heat 
data. In what follows we first present the main theoretical 
framework, and then report a detailed analysis of the experi- 
mental data within this framework. The central result of our 
analysis is the doping dependence of the pairing pseudogap 
energy extracted from specific heat data, which compares well 
with the doping dependence of the pseudogap scale inferred 
from a variety of other measurements. 

Pairing contribution. For the calculation of the pairing 
contribution we start with the usual expression of the grand 
potential for pair fluctuations [ |l4| ] 

n p (T) = —T ln [1 - 9X (q, w»)] , (1) 
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where x(q, ui n ) is the pairing susceptibility. In what follows 
we will consider only the non-interacting, bare susceptibility: 



x(q, u n 







+ q 1 iCm+iu>n)GP(k,iCm): ( 2 ) 



where [G°(k, «Cn)] _1 — Kn — £k- Here £k is the single parti- 
cle dispersion, £ n = (2n + l)7r//3 is the fermionic Matsubara 
frequency. The summation in Eq. ([]]) is done over the pair 
momentum q, bosonic Matsubara frequencies oj n = 27rn//3 
and spin a. We consider the renormalized classical regime 
corresponding to uj„ ~ 0. Classical fluctuations give in two- 
dimensions the dominant contribution to the self-energy at 
low frequencies. Quantum fluctuations (ui n ^ 0) are impor- 
tant only at low temperatures and are irrelevant for temper- 
atures T c < T < T*. Thus, in the renormalized classical 
regime the susceptibility is given by the relation: 



2„2 



1 - 9X (q) = e T + £o4 



(3) 



The first term on the right side of the above expression 
corresponds to an exponentially growing correlation length, 
£ ~ exp(— T*/T), a hallmark of the classical fluctuation 
regime. The onset of this regime happens at T*, a charac- 
teristic temperature scale set by the coupling strength g and 
pseudogap A (note that in this notations g has units of en- 
ergy): 



k B T* = —A 2 . 
9 



(4) 



For the present purposes we identify the self-energy S(k, ui n ) 
in the above expression with the pairing self energy JTo| ] of 
the renormalized classical regime: 



S(k,w n ) 



A 2 



(7) 



where A is the pseudogap energy scale. For the dimension- 
less specific heat coefficient coming from this contribution we 
obtained the result: 



7-(T) 



2tt 2 



N(x,d) 
N(0) 



1 ) f(x)(l-f(x))x 3 



(8) 



T* can also be identified with the onset of pseudogap 
behavior. In principle both A and g can be doping 
dependent. The dimensionless specific heat coefficient 
7(T) = c v (T)/k 2 B TN(0) can be obtained directly from 
Eq.(l) according to the thermodynamic relation c v (T) = 
—Td 2 ft I dT 2 . The result for the pair contribution is 

1p {T) = X^a{^\ e- 2 ^ fl + JLi nA Y (5) 

In strictly two dimensions the integration of the pair mo- 
menta needs to be regularized by introducing an upper cutoff 
A. In practice, however, the materials we are interested in 
have a large, but finite, c-axis anisotropy, that naturally cuts 
off the pair momenta A ~ £o /d c where d c is the interlayer 
distance. In Eq. (§) A ~ {Qk B T*N (0)) _1 is the parameter 
that governs the relative importance of the pairing contribu- 
tion. The characteristic feature of the pairing contribution to 
the specific heat is a broad hump near a temperature T ~ T* . 
This contribution in general, and the hump feature in particu- 
lar, are important for obtaining a good fit to the experimental 
data. 

Single particle contribution. The single particle contribu- 
tion to the grand potential is given by 

il s (T) = -T { ln t G_1 ( k > w «)] -£(k,0 • G{k,u n )} 



In Eq.(||) we introduced the following notations: x — 
E*/k B T, d = A/k B T, E k = ^/ei + A 2 , f (x) = 
(e x + 1) . As indicated by angle-resolved photoemission 
spectra [|3|] taken in the pseudogap regime, the single particle 
states in this regime are very broad. During the analysis of the 
specific heat data we found that even for near-optimal doping 
an intrinsic broadening of the single particle states must be 
considered in order to obtain a reasonable account of the data 
based on the single particle contribution (||). The importance 
of broadening in the single particle states becomes even more 
apparent when our analysis is performed on data taken on un- 
derdoped samples. This trend is again consistent with that 
observed in ARPES. We therefore introduce a phenomeno- 
logical broadening of the single particle density of states in 
the following simple form 



N(x,d) 
N(0) 



Re 



^7 



J (x — irff + d 2 



(9) 



and 7 = T/k B T (not to be confused with j p (T) or 7 S (T)). 
We would like to note that while the broadening V is essential 
for treating the single particle contribution it is less important 
for the qualitative behavior of the pair contribution 7 P (T) as 
function of temperature. The total contribution to the specific 
heat coefficient is 



lel {T)= ls {T)+ lp {T). 



(10) 



(6) 



Relative magnitude of the pairing and single particle con- 
tributions. From the results obtained in Eqs. (4) and (7) 
we see that the relative importance of pairing contribution 
7 P (T) compared to the usual contribution 7 S (T) is governed 
by the parameter A - {^k B T* N . Using £ ~ 10 A 
T* ~ 100K and typical density of states, (for example, from 
ARPES data we get 3 states /eV per Cu — O in a unit cell 
[j^]), we find that j p and 7 S are comparable in magnitude: 
7 P (T) ~ Js(T). In contrast, for a superconductor with long 
coherence length such as Al (^ai ~ 100£ybco) the pairing 
contribution j p (T) is entirely negligible. 

Comparison with the experiment. A typical comparison of 
the experimental data with the theoretical result is presented 
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in Fig.l. The single particle and the pairing contribution is 
shown separately in order to emphasize the relative size and 
the importance of both terms. 




300 



FIG. 1. Fit to the experimentally measured specific heat of 
Y BiCuzOn.iz (dots) using the theoretically calculated single par- 
ticle and pairing contributions (lines). Note that the single particle 
contribution "/ s (T) cannot reproduce the hump and the pairing term 
7 P (T) is necessary. 



stant is obtained as described above, [ 0.14 mJ/gatK 2 for 
La2-xSr x CuOi and 0.23 mJ/gatK 2 for Y Ba^Cu^O^x 
] it is no longer a fitting parameter and remains fixed for 
all doping levels. When considering both single particle and 
pair contributions, a direct comparison of the theoretical re- 
sults and the experimental data can be performed. The results 
are presented in Fig.2 for Y Ba^Cu^O^x, and in Fig.3 for 
La 2 -xSr x Cu0 4: 0. 
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The single particle contribution monotonically increases 
with increasing temperature and by itself cannot explain the 
hump observed in the experimental data. If, however, the pair- 
ing contribution j p (T) (with maximum near T*) is included, 
the fit is reasonable [O]. During the fitting procedure only the 
normal state (T > T c ) specific heat values are used. 
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FIG. 2. Specific heat data (dots) for Y B-zCu^Os+x for a wide 
range of doping (from below) x = 0.43, 0.48, 0.57, 0.67,0.73,0.76 
and 0.80, and the corresponding theoretical curves (lines). The inset 
gives extracted values for A and V. 



FIG. 3. Same as in Fig. 2, but now for La,2-xSr x CuOi 
[x = 0.03, 0.05, 0.084, 0.1, 0.125, 0.135 and 0.15]. 

The phase diagram obtained from the extracted values of 
the pairing pseudogap A, and the experimental transition tem- 
peratures, is presented in Fig A. Here we use reduced quan- 
tities (normalized to their value at optimal doping) to present 
on the same scale the experimentally measured critical tem- 
peratures T c of Y Ba^CuzOs+x and La2- x Sr x CuOi com- 
pounds, and the corresponding pairing pseudogap scale A ex- 
tracted from data. 
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An overall factor of order of 0.1 and dimensions of 
mJ/gatK 2 is used to scale the dimensionless theoretical 
expressions to the data at optimal doping. Once this con- 



FIG. 4. Doping dependence of the pairing pseudogap A ex- 
tracted from specific heat and of the critical temperature T c for 

YBa2Cu30e+ x and La2- x Sr x Cu04,. 
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An interesting result emerges from the present analysis. 
The best fit is obtained if T*(x) is taken to be proportional 
to A(x): T* ~ A(a;). 
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FIG. 5. T* and T c dopind dependence measured as 

ARPES, 7, x for YBa 2 Cu 3 6 + x , La 2 - x Sr x CuOi and 
B12 Sr 2 CaCu-i Og+ x . 

This proportionality has also been indicated in Ref. [ |l7| ] 
based on a scaling analysis of various experimental data, and 
an approximate ratio of A = 1.5T* has been found, similar to 
our finding. Eq.(3) implies, at least for the present pair fluc- 
tuation formalism, that the ratio A(x) /g{x) is independent of 
doping. Thus the doping dependence of the A(x) / A(x max ) 
coincides with doping dependence of T* {x) / T* (x max ) . The 
broadening T(x) increases with decreasing the doping, which 
is again in qualitative agreement with the ARPES results [|^]. 

Conclusions. We presented the first systematic theoretical 
analysis of the electronic specific heat in the pseudogap state. 
While this analysis was performed within the classical pair 
fluctuation framework for the pseudogap, we suspect that the 
general features of the results are relevant for a large class of 
pairing fluctuation scenarios. 

The main results of this analysis can be summarized as fol- 
lows: (i) both the single particle and pair contributions are 
of the same order of magnitude and needed to fit the specific 
heat data of pseudogapped samples, (ii) the electronic spe- 
cific heat data of underdoped cuprate samples have several 
universal features that are well captured by our pair fluctu- 
ation model. Our analysis of the experimental data yielded 
a pairing pseudogap energy scale A(x) that has a magnitude 
and doping dependence which is similar to the pseudogap en- 
ergy scales extracted from other experiments (cf. Fig.5). Our 
main result is the phase diagram presented in Fig A. 
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